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This article provides a gentle introduction for a general mathematical audience to the factorization theory of motion 
polynomials and its application in mechanism science. This theory connects in a rather unexpected way a seemingly 
abstract mathematical topic, the non-unique factorization of certain polynomials over the ring of dual quaternions, with 
engineering applications. Four years after its introduction ©[TOl, it is already clear how beneficial it has been to both 
fields ||6 l[T2l (17112^ . In Section [T] we infroduce the notion of motion polynomials and discuss their decomposition 
into products of linear motion polynomials. It can be used to synthesize linkages following a prescribed motion and is 
related to a variant of Kempe’s Universality Theorem. We explain the relation to mechanism science in more detail in 
Section]^ In Sections[^and[^we present examples from linkage synthesis and discuss exceptional factorizations. 

1 Motion polynomials and their factorizations 

We denote by El the skew field of quafemions, generafed by i,j,k over M, with the well-known relations 

P =j2 = k2 =ijk= -1. 

The conjugate of a quaternion q = q\+ q 2 \ + q^\ -|- <74k is defined by ^ = <71 — q 2 i — <?37 — q^k, the norm of q is 
N{q) := qq = q\ + q\ + q\ + ql- The scalar extension DEI := D (g)® El by dual numbers D := M[e]/ {e^) is just a skew 
ring, the ring of dual quaternions. The conjugate of a dual quaternion h = p + eqi^ defined hy h = p + Eq. We also 
use N{h) to denote the norm of h with N{h) := hh = N{p) -f £{pq + qp) G D. The dual quaternion h is invertible if 
and only if p 7^ 0. 

Denote by S the multiplicative subgroup of dual quaternions with nonzero real norm. Its elements may be written as 
h = p + £q where the quaternions p and q satisfy p ^0, pq + qp = 0. The latter equality is called the Study condition. 
The group S acts on = (i,j,k) according to 


X I—)• 


pxp + pq-qp 

W) 


( 1 ) 


Any such map is a Euclidean displacement. The rotational part is the well-known action x 1—)• pxp/N{p) of the unit 
quaternion pA(p)^72 on M^, the translational component is {pq — qp)/N{p). Equation o defines a homomorphism 
from S fo fhe group SE(3) of Euclidean displacemenfs. If is surjecfive and ifs kernel is fhe real mulfiplicafive group M*. 

This algebraic construction allows a geometric interpretation. Identify DEI/M* with real projective space of 
dimension seven, denote by ^5^ C : pq + qp = 0 the Study quadric and by E the three-space with equation p = 0. 
Then Study’s kinematic map is the bijection SE(3) = S/E* 5^\E C whose inverse maps h = p + eq to the 
Euclidean displacement given by O- 
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The concept of motion polynomials arises by making the above group homomorphism parametric. More precisely, 
let D]HI[f] be the skew ring of univariate polynomials over DEI in one variable t that commutes with all coefficients. For 
C G D]HI[f], the conjugate polynomial C is obtained by conjugating all coefficients. If C = P + eQ with P,Q^ then 
we call P the primal part and Q the dual part. We say that C G D]HI[t] is a motion polynomial if N{C) := CC G M[f] \ {0} 
(a priori, the norm polynomial N{C) is in D[t]) and if its leading coefficient is invertible. For any to gM. which is 
not a zero of N{C), we can say that C(to) is an element of S, hence it acts on M^. Varying to, we get a motion, i.e. a 
parametrized curve of Euclidean displacements. 

By virtue of ([T]), the trajectories of points during that motion are rational curves whence the motion itself is called 
rational. It is well-known that all motions with only rational trajectories have a polynomial parametrization with 
values in the Study quadric lfT4l . The motion polynomials do not constitute a multiplicative group. However, the 
quotient by M[f] — { 0 } is a group because the inverse of the class of a motion polynomial is precisely its conjugate. 
While kinematic or geometric properties of motion polynomials do not change if we multiply them with non-zero real 
polynomials, algebraic properties may be different. This observation will be important in Section]^ 

Summarizing, we can state 

Proposition 1. Motion polynomials parametrize rational motions and the group of motion polynomials modulo 
M[t] \ {0} is isomorphic to the group of rational motions. Via Study’s kinematic mapping, motion polynomials (rational 
motions) correspond to rational curves on the Study quadric SP with at most finitely many points in the three space E. 

Since we want to use some version of the fundamental theorem of algebra later on, it is convenient to restrict 
attention to monic motion polynomials. As far as applications in kinematics are concerned, this is no loss of generality 
and can always be accomplished by a suitable coordinate change. The fundamental theorem speaks about factorization 
into linear polynomials, so it is time to clarify the kinematic nature of linear motion polynomials. 

Proposition 2. Every monic linear motion polynomial parametrizes either a rotation about a fixed axis in or a 
translation in a fixed direction in M^. 

The converse is not true: there are monic motion polynomials parametrizing a rotation around a fixed axis that are 
not linear. 

Example 1. The dual quaternion polynomial Ci =t — i has norm t^ + 1, hence it is a motion polynomial. It parametrizes 
a rotation around the first coordinate axis. Indeed, by ([T]l we have 


(t^ — \ 2t \. —I 2t \, 

Setting (p := 2arctant, i.e., costp = — t^)/(t^ + \), sintp = 2t/(t^ + \), the assertion becomes apparent. 

The dual quaternion polynomial C 2 = (f — i)^ = — 1 — 2it has norm (t^ + 1)^, hence it is a motion polynomial. It 

also parametrizes a rotation around the first coordinate axis. Indeed, it can be obtained by reparametrizing the first 
parametrization using the reparametrization function 11 —)• — I)/(2t) and clearing denominators. 

Example 2. The dual quaternion polynomial C 3 = f — £i has norm t^, hence it is a motion polynomial. The action on 
is 

2 ei 

xii-I-X 2 J -|-X3k I—>■ xii-\-X2j -|-X3k H ——. 

This is a translation in the direction parallel to the first coordinate axis. 

The following theorem is fundamental in the factorization theory of quaternion polynomials |7]: 

Proposition 3. Every monic polynomial C G admits a factorization C = (t — qi) ■■■ (t — qn) with quaternions 
^1 ) ■ • • ) (in- 


X\\+X-2^+X-i\i\ 


(t-i)(xii-hX2j+.r3k)(t-ki) 
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For the proof, one uses a fundamental theorem for monic non-negative real polynomials: any such polynomial (it 
must he of even degree) has a unique factorization into monic non-negative real quadratic polynomials. This is an easy 
consequence of the fundamental theorem for complex polynomials. The norm polynomial A/^(C) is a non-negative 
real polynomial and factors into the product of non-negative real quadratic polynomials. The Euclidean algorithm for 
polynomial division also works in the case of polynomials in EI[t]. If we take one of the non-negative quadratic factors 
of the norm, say 2 G ]R[f], then the polynomial remainder of C hy 2 is a linear quaternion polynomial or zero. In the 
first case, this linear polynomial can he factored out. In the second case, Q can he factored out. 

A specialty in this context is that a generic polynomial C G EI[f] has precisely n\ different factorizations. Since 
the ring of quaternion polynomials is not commutative, permuting the factors no longer yields a factorization C. In 
other words the n ! factorizations differ essentially, whereas in the complex case all factorizations can he obtained as 
permutations of one. In this case we rather would think of just one factorization whose factors can he permuted. 

Another word of caution concerns the polynomial division mentioned in the above sketch of the proof. Since ]HI[f] is 
not commutative, we have left and right polynomial division, left and right quotients, and left and right remainders. It 
does not matter which direction we choose, but we have to choose one and stick to it. In the following we deal only 
with right polynomial division, right remainders, and consequently left quotients: (C = Quotient • D + Remainder). 

Let us pass from polynomials over ]HI[f] to motion polynomials in DEI[f] and try to factor the motion polynomial C. 
Like in the proof above, the norm is non-negative and, by the important defining property of motion polynomials, 
real. The quotient by one of the quadratic factors has degree at most one. But now, polynomial division might be 
problematic because the linear remainder polynomial is in general not monic and the leading coefficient might not 
be invertible. It is even possible that the remainder is constant, as in the example C = t^ + \ + e\ where we have 
CC = + 1)^ and the remainder is ei. 

However, in the generic case, everything is fine: the remainders will be linear with invertible leading coefficients 
and can be divided out. We can then construct step by step a factorization into linear factors by Algorithm [T] below. 
This genericity condition can be simplified as follows: For Algorithm [T] to work, it is sufficient that the primal part of 
C has no nontrivial real factors. The different factorizations come from the arbitrary choice of a quadratic factor M in 
Line 1^ of the algorithm. 


Algorithm 1 (factorization of generic motion polynomials) 

Input: Motion polynomial C, monic, primal part has no real factor. 

Output: List L= [fii,..., /i„] of generic motion polynomials 



such that C = {t — hi) ■■■ (t — hf). 


1 

F ^ list of quadratic, irreducible factors of CC 


2 

D^C 


3 

initialize L = [ ] 

0 empty list 

4 

while F is not empty do 


5 

choose M £ F 


6 

write D = QM + R with degR < 1 

0 polynomial right division 

7 

compute unique zero hof R 


8 

prepend htoL 


9 

D ^ D' where D = D'{t — h) 

0 polynomial right division 

to 

delete M from F 


11 

end while 



In fact, we have 

Theorem 1 ( IIIOll ). Algorithm^can be used to factor a motion polynomial C = P + sQ, provided the primal part P 
has no real factors. In this case, C = {t — hi) ■■■ f — hf) and for i G {I,... ,n} each polynomial t — hi describes a 
rotation about a fixed axes. Moreover, all possible factorizations (in general n\) can be obtained in that way. 

As a consequence of Theorem [T] we have 
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Figure 1: Anti-parallelogram linkage 


Corollary 1. In general, a rational motion of degree n can be decomposed in at most n\ different ways into the product 
of rotations t — h\, ..., t — h„. 

The phrase “in general” in Corollary [T]refers to the absence of real factors in the primal part of C. If this requirement 
is not fulfilled, the statement is not true. In fact, all kinds of special behavior can be observed. It is interesting that these 
“exceptional” examples surprisingly often arise in natural applications of Theorem [T] to kinematics and mechanism 
science. We will return to this point a little later. But at first, we have to explain the relation between Theorem [T] and 
mechanism science. 

2 Factorizations and linkages 

Consider a generic motion polynomial C of degree two that parametrizes a planar motion (all trajectories are in the 
parallel planes). By Theorem[T] it admits two factorizations 


C = {t — h\){t — hi) = {t — k\){t — ki) 

with suitable dual quaternions hy^hi^ki^ki. This means that the motion parametrized by C can be generated in two 
ways as composition of two rotations with respective centers h\,hior k\,k 2 . Moreover, the two revolute joints at hi 
and ki can be rigidly connected without disturbing the motion C. This is illustrated in Figure which depicts a planar 
four-bar linkage whose joints are labeled by the corresponding dual quaternions. The joints hi and k\ are fixed, hi 
and ki can rotate about hi and ki, respectively, but retain their distance. Planar four-bar linkages constitute the most 
important class of linkages for engineering applications. Our case is special because the joints hi, hi, ki, and ki form 
an anti-parallelogram. 

If the motion polynomial C is not planar, a similar construction yields a spatial linkage, consisting of a closed 
loop of four skew revolute axes (Figure]^. Any two neighboring axes are rigidly connected, that is, they maintain 
their distance and angle throughout the motion. In contrast to the planar case, the one-dimensional mobility of such 
a structure is not obvious. A closed loop of four revolute joints is generically rigid. But the loops resulting from 
factorizations of quadratic motion polynomials move because of their algebraic construction. 

Spatial four-bar linkages that move with one degree of freedom have been known for a long time 121 12 [161. There 
exists only one family of this linkage type and its members are referred to as “Bennett linkages”. So far, Bennett 
linkage have minor importance in applications but we shall see their theoretical significance later in this text. Moreover, 
they exhibit a fascinating geometry. A stunning example is the Wunderlich’s explanation of the Backlund transform of 
discrete asymptotic nets of constant Gaussian curvature by means of Bennett linkages l28l . 

By factorizing cubic motion polynomials, we can also construct spatial six-bar linkages with a one-dimensional 
mobility as in Figure Their complete classification is a long-standing open problem in theoretical mechanism 
science. In spite of some recent progress |[8l [TTl, a classification is currently still out of reach. At any rate, our 
approach yields new examples of spatial six-bar linkages IITtI - IT^ and the first viable synthesis procedure. We describe 
this in more details in the next section. 
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Figure 2: Three-pose synthesis of a four-bar linkage 


3 Linkage synthesis 

An important application of motion polynomial factorization is linkage synthesis, the construction of a mechanical 
linkage to accommodate a certain task. Our approach is well suited to exact synthesis with prescribed poses, that is, 
the computation of linkages such that one link visits a finite number of prescribed poses. The word “pose” refers here 
to the position and orientation of a rigid body, that is, an element of SE(3). 

Let us consider a simple example of a spatial four-bar linkage (Bennett linkage). Its coupler motion is given by 
a quadratic motion polynomial C. In the kinematic image space , it parametrizes a conic section on the Study 
quadric Conversely, a generic conic section on y gives rise to generic quadratic motion polynomials (differing 
only by admissible reparametrizations) and, via motion factorization, to a Bennett linkage. Thus, we can synthesize a 
Bennett linkage to three prescribed poses po, p\, p 2 - These poses are points on the Study quadric y where they span 
a plane. We compute a rational quadratic parametrization C of the intersection conic of this plane and y and factor it 
as C = (t — h\){t — h 2 ) = (t — k]){t — ^ 2 ). The linear motion polynomials in these factorizations determine the fixed 
axes {h\, ki) and the moving axes {h 2 , ^ 2 ), as illustrated in Figure]^ 

Our synthesis procedure for Bennett linkages is elegant and simple but it is not the only available method. Other 
approaches include ||4 tl24l - l27ll . It extends, however, to linkages with more than four joints. We may, for example, 
synthesize six-bar linkages to four prescribed poses |[T^ . Here, the geometry is slightly more involved but still 
understandable via classical results. Four points po,pi,p 2 ,P 3 in general position on the Study quadric y span a 
three-dimensional projective space. If this space intersects ,5^ in a ruled quadric, as in Figure]^ the four points can be 
interpolated by two one-parametric families of rational cubic curves. Each interpolating cubic can be factored and 
gives rise to several open chains with three revolute joints that can be combined to form a closed six-bar linkage. Here 
are a few further remarks on this construction: 

• Closed loops of six revolute axes are, in general, rigid. In our case, they move because of their special algebraic 
construction. 

• There exist spatial six-bar linkages that have a one-parametric mobility but cannot be constructed by factorization 
of motion polynomials. The relative motions between two links are not all rational or, at least, do not have 
rational components. 

• The above construction is the first viable synthesis procedure for spatial six-bar linkages. It may easily be 
generalized using a Hermite like interpolation scheme. 


4 Exceptional factorizations and Kempe’s Universality Theorem 

Theorem [T] shows the existence of finitely many factorizations of generic motion polynomials. In this section we are 
concerned with non-generic situations where the primal part of the motion polynomial has real factors. In this case. 
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Figure 4: Three different factorizations of a circular translation 


Theorem [T] gives no information. The following examples demonstrate what can happen: 

Example^. The motion polynomial C = (t — l)(t — j) — £((i + k)t — 2k) can be factored as 

C = (f — 1 — £i)(t — j — £k) = (t — j — £(i + 2k))(t — 1 +£k). 

The polynomial factors f — 1 — £i and f — 1 — £k parametrize, however, translations, not rotations. One may view this 
as a limiting case of the generic situation appearing in Theorem [T] in the sense that the two rotations degenerate to 
translations. It turns out that they can still be computed by Algorithm[T] 

Example 4. We consider the motion polynomial C = + 1 + £(ai + b\t) with a^b a^b >0, + b^ > 0. It 

parametrizes the curvilinear translation along an ellipse with semi-axis lengths a and b. \i b, a straightforward 

computation shows that C admits no factorization of the form C = {t — h\){t — h 2 ) with linear motion polynomials 
t — h\ and t — /i 2 . If, however, a = b (circular translation), even infinitely many factorizations exist, namely 

/ji =k-£(/i+(a + g)j), /i2 =-k + £(/i + gj), (2) 

with /, g G M. They have a very clear geometric explanation. The motion in question is a circular translation and can 
be generated by infinitely many parallelogram linkages (Figure |^. Each leg corresponds to one factorization of the 
form Q. 

Existence of exceptional situations demonstrate that the factorization theory over DB[[t] is more complicated but 
also more interesting than the theory over EI[f]. Moreover, situations with no or infinitely many factorization arise 
surprisingly often in engineering applications of motion polynomial factorization. Many important rational motions 
are not amenable to straightforward factorization via Algorithm [T] Still, it is possible to factor these motions but at the 
cost of raising the number of factors (the degree of the motion polynomial). Recall that for any non-zero polynomial 
/? G M[f] the motion polynomials C and CR parametrize the same motion. Thus, one may try to find a real polynomial 
R such fhaf CR admifs a facforizafion. 
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Figure 5: Linkage to generate an elliptic translation and corresponding link graph 


In engineering applications, revolute joints are preferred over translational joints. Hence, we focus on factorizations 
with revolute joints only. In this case, a necessary requirement is that the motion parametrized hy C is hounded, i.e., all 
trajectories are bounded rational curves. We also say that the motion polynomial C itself is bounded. Indeed, we have 

Theorem 2. For every bounded motion polynomial C there exists a real polynomial R such that CR admits a 
factorization CR = {t — hi) ■ ■ ■ {t — h^) with rotation polynomials t — hi, t —hn- The degree ofR is bounded by the 
maximal degree of a real factor of the primal part ofC. 

Theorem|^has been proved in 10 for the planar case and in ll^ for the general case. The proofs are constructive 
so that factorizations can be effectively computed. The main ingredients are variants of the Euclidean algorithm for 
polynomial division over the (dual) quaternions and the solution of quadratic equations with real coefficients over 
quaternions lfT3l . It is worth noting that infinitely many factorizations exist if deg/? > 0. 

As an application of Theorem]^ consider the elliptic translation appearing in Example]^ There exists a quadratic 
real polynomial R such that CR admits the factorization CR = {t — hi){t — h 2 ){t — hf}{t — hf). From this, we construct 
an (admittedly complicated) linkage to generate an elliptic translation (Figure [^: 

• The linkage consists of four anti-parallelograms for i G {1,2,3,4}). 

• Six angles (<(/i,,m,',/j,+i), for i G {1,2,3}) are kept constant during the motion. In other words, 

we have a chain of anti-parallelograms, where each anti-parallelogram “follows” its predecessor. 

• The rigid body attached to the connection of m 4 and h 4 performs the elliptic translation. The point h 4 “draws” 
the indicated ellipses. 

Figure also shows a more abstract representation of the same linkage. In this “link graph”, each vertex represents a 
link (a rigid connection between revolute joints) and each edge represents a joint. Two vertices are connected, if the 
corresponding links share a joint. The linkage of Figurej^was constructed by augmenting the given factorization (hi, 
h 2 , /i 3 , ( 14 ) with one additional joint mo and then successively computing the remaining joints by solving the recursion 

{t-mi^i){t-hi) = {t-ki){t-mi), /G {1,2,3,4} (3) 

with the aid of Algorithm[T] We call this computation of ki and m, from and a Bennett flip. This name comes 
from the fact that in the spatial case the involved dual quaternions determine the axes of a Bennett linkage. In the 
planar case, a Bennett flip generates anti-parallelograms. 


7 






Figure 6 : Linkage to draw an ellipse and link graph 


Figure [^presents a refinement of this construction. Instead of multiplying C with a polynomial G M[f], we right 
multiply it with the linear polynomial // := t — k G IH[[t] such that the product CH admits a factorization. Of course, 
C and CH parametrize different motions but the fixed points of H (points of the third coordinate axis) retain their 
trajectories. Flence, we obtain a linkage that is capable of drawing the prescribed ellipse. The triples of collinear 
joints come from the fact that the constant angles in the linkage of Figurej^are straight for the linkage of Figure]^ 
The motion of the link connecting m 3 and /13 is no longer an elliptic translation and generic trajectories are of degree 
greater than two. 

It is a yet unpublished result that one can always find a polynomial H G such that CH admits a factorization. 
Using the Bennett flip technique, one can then construct a (in general) spatial linkage with prescribed trajectories. 
The anti-parallelograms of Figure [^become Bennett linkages and the additional freedom we gain by multiplying 
with a quaternion polynomial // G IHI[t] instead of a real polynomial G M[f] can be used to reduce the number of 
anti-parallelograms. 

The constructions of this section are not restricted to ellipses but can be generalized to arbitrary rational space 
curves. First let D denote a rational curve of degree d in the Euclidean 3-space. Then consider a motion polynomial 
C that parametrizes a rational motion with trajectory D, for example the curvilinear translation along D. Using the 
algorithm of Theorem]^ we arrive at a factorization of CR as {t — h\).. .{t — h,„). This gives us an open chain of m 
revolute joints representing our motion polynomial C. Using Bennett flips as in Q, we can construct linkages with 
only revolute joints that draw arbitrary (bounded) rational curves. 

This is not a new result. By a celebrated theorem of mechanism science (Kempe’s Universality Theorem ||5j[T5]ll. 
any bounded portion of an algebraic curve can be drawn by a linkage. Our approach shows that in the rational case the 
number of necessary links and joints decreases dramatically. The asymptotic bound for rational curves is linear in the 
curve degree n while the currently best known bound for general algebraic space curves is cubic l|T]. Moreover, our 
general construction behaves quite well in important low-degree cases. The linkage in Figure |^has only ten joints 
while the ellipse linkage in Kempe’s construction requires as much as 235 joints. Often it is possible to further reduce 
the number of joints so that engineering applications come into reach. For example, only seven joints are necessary 
to generate the so-called Darboux motion, a spatial motion where all trajectories are ellipses in non-parallel planes 

inniEii. 
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5 Conclusion 


Motion polynomial factorization is an algebraic theory with surprising relations to mechanisms science. The interplay 
between both disciplines is beneficial to both. Algebra can provide solutions to hitherto inaccessible engineering 
problems and requirements of applications led to interesting algebraic questions. Our current work focuses on both, 
the details of a version of Kempe’s Universality Theorem for rational space curves with emphasis on a low number of 
links and joints and further applications of motion polynomial factorization to engineering problems. 
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